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STRESS  CONCENTRATION  ABOUT  CURVILINEAR  HOLES  IN 
PHYSICALLY  NONLINEAR  ELASTIC  PLATES 


**71009 


A.N.Guz1,  G.N. Savin,  I.A.Tsurpal  (Kiev) 


An  approximate  solution  method  of  plane  physical  nonlinear 
problems  of  ..stress  'concentration  about  curvilinear  holes 
in  thin  plaxes  rafle  of  a  material  subject  to  a  nonlinear 
law  of  elasticity  is  given.  The  solutions  are  represented 
in  the  form  of  expansions  in  the  small  parameter  |i  and  e. 

The  determination  of  the  stress  function  F  for  a  physically 
nonlinear  elastic  plate  with  a  hole  reduces  for  each  ap¬ 
proximation  to  the  integration  of  nonlinear  differential 
equations.  Stress  concentration  about  an  elliptic  hole  is 
considered  in  zero,  first,  and  second  approximation.  The 
coefficient  of  stress  concentration  k  is  found/pn  the  contour 
of  the  hole,  depending  nonlinearly  on  the  tensile  forces  P, 
the ellipticity  of  the  hole,  and  a  parameter  X  characterizing 
the  mechanical  properties  of  the  material.  Tables  represent 
the  values  of  the  coefficient  of  stress  concentration  for 
various  values  of  the  parameters  P,  X ,  and  e .  I  . 


1.  The  problem  considered  is  that  of  stress  concentrations  in  the  neighbor¬ 
hood  of  curvilinear  orifices  without  sharp  corners  in  a  thin  plate  consisting 
of  a  material  for  which  the  stress  -  strain  ratio  is  nonlinear  even  in  the 
presence  of  comparatively  small  strains.  Given  the  deformation  magnitude  con¬ 
sidered,  all  geometric  relations  of  elasticity  remain  linear,  i.e.,  we  are  deal¬ 
ing  with  a  variant  of  the  physically  nonlinear  theory  of  elasticity,  with  a 
specified  nonlinear  law  of  elasticity. 

A  previous  study  by  the  author  (Bibl.l)  examined  this  problem  for  the  non¬ 
linear  law  of  elasticity*  (Bibl.2),  using  conformal  napping  of  the  region  in 
question,  outside  the  curvilinear  orifice,  onto  the  exterior  of  a  unit  circle 
and  introducing  the  Kolosov-Muskhelishvili  complex  potentials.  For  the  sought 
stress  function,  represented  in  the  form  of  expansions  in  a  small  parameter, 
differential  equations  and  boundary  conditions  for  successive  approximations  in 
curvilinear  coordinates,  given  by  a  mapping  function,  have  been  derived.  How¬ 
ever,  in  view  of  the  embersomeness  of  the  right-hand  sides  of  the  equations, 
this  method  has  led  to  extremely  complex  calculations  with  respect  to  orifices 
of  noncircular  shape. 


This  paper,  utilizing  the  same  nonlinear  law  of  elasticity  (Bibl.2),  pro- 


*  The  stress  concentration  in  the  neighborhood  of  a  circular  orifice  for  this 
law  of  elasticity  has  been  investigated  elsewhere  (Bibl.2,  7,  8,  and  9). 

**•  Numbers  in  the  margin  indicate  pagination  in  the  original  foreign  text. 


poses  another  approximate  method  for  the  solution  of  the  above  problem,  which 
makes  it  possible  to  complete  this  solution  with  respect  to  certain  noncircular 
orifices.  The  new  approach  is  based  on  the  approximate  method  of  "perturbation 
of  the  boundary  form"  (Bibl.3),  as  successfully  used  by  the  authors  (Bibl.A,  5, 

6)  in  investigating  stress  concentrations  in  the  neighborhood  of  analogous  curvi¬ 
linear  orifices  in  shells. 

2.  The  approximate  method  described  here  requires  the  representation  of  all 
basic  equations  in  a  polar  coordinate  system;  hence  we  will  use,  in  the  form 
given  elsewhere  (Bibl.2),  the  nonlinear  law  of  elasticity  for  a  generalized 
plane  stress  state; 


ef  =  -—*(«,)  or„  +  ~gOS)((Tr  “  <To) *  /1010 

^=‘jj^k(s0)a0  +  -~-g(tS)(<rr-(ro).  (2.l) 

erp  —  ~q~  8  (tl)  » 

where  er  ,  e^,  and  ercp,  as  well  as  crr  ,  csy,  and  TrCp,  correspondingly,  are  the  mean 
stress  and  strain  components  over  the  plate  thickness  in  the  polar  coordinate 
system  (r,  cp) ;  K  and  G  are  constant  moduli  of  volume  deformation  and  shear, 
respectively,  for  the  physically  nonlinear  material  of  the  plate  in  the  presence 
of  vanishingly  small  deformations;  k(so)  and  g(t^)  are  the  pressure  and  shear- 
stress  functions  which  characterize,  respectively,  the  change  in  volume  and 
shape  at  any  point  of  the  body  during  its  deformation.  The  dimensionless  quanti¬ 
ties  s0  and  to  are  expressed  in  the  form  of  invariants; 


S°  3K  “  9K  + 

t3  =  =  <4*^* +a*~  a,°* + 3r*^  • 


(2.2) 


For  many  materials,  volume  deformation  over  a  wide  range  obeys  Hooke *s  law  so 
that,  in  eq.(2.l),  it  may  be  assumed  with  a  high  degree  of  accuracy  that 
k(3o)  -  1.  The  slight  deviations  from  linear  dependence  between  stress  and 
strain  in  the  elasticity  relations  (2.1)  can  be,  with  sufficient  accuracy, 
mapped  by  the  function  g(t^)  =  1  +  ggt^  . 

Hereafter,  we  will  assume  that,  in  the  elasticity  relations  (2.1), 


*(*.)  —  1  >  S(<o) — 


(2.3) 


where  &  is  a  dimensionless  constant. 

At  such  a  choice  of  the  nonlinear  law  of  elasticity  (2.1)  under  the  condi¬ 
tions  (2.3),  the  problem  of  the  stressed  state  of  a  thin  plate  reduces  to  find- 
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ing  the  stress  function  F(r,  cp)  from  the  fourth-order  nonlinear  differential 
equation  (Bibl.2,  8)  : 


AAF+  lf\yF'T*  +  F*TrJ 

-  -^-(-jFrrT,*  -  Fr,T„  +  yF„r„) -  ~{FtT„  +  F„Tt)  -  y A  {. TdF )]  =  a 


(2.4) 


in  the  presence  of  corresponding  boundary  conditions  over  the  orifice  contour 
and  at  ’•infinity”.  If  the  function  F(r,  cp)  is  known,  the  stress  components  oT  , 
crcp,  and  TrCp  may  be  directly  found  with  the  aid  of  F(r,  cp)  from  the  formulas 


1  r  1  dF  ,  1  dlF~L  1  r  1  r*  ,  1  „  1. 

a'  -  R*  l  r  dr  +  r *  dtp*  J  “  R*  [  r  F' +  r*  F"]1 


(2.5) 


In  eqs.(2.4)  and  (2.5)»  r  is  a  dimensionless  coordinate  referring  to  the  /1011 
quantity  R  which  characterizes  the  absolute  dimensions  of  the  orifice;  A  is  the 
Laplace  operator  which,  in  dimensionless  coordinates,  has  the  form 


A  =  ~Fr  +  -prFff+F„. 

The  material  constant  X  and  the  function  T(r,  cp)  have  the  form 


(2.6) 


Kgt 


(3K+G)Gl 


(2.7) 


where  u  =  — - —  is  a  small  dimensionless  quantity 


«.  =  K*s  - 

P  (3K  +  G)Ga’ 


T(r,<p)  —  yGaf* » 


where  t^  is  given  by  the  expression  (2.2)  and  the  components  crr  ,  a„,  and  Tfy  are 
associated  with  the  stress  function  F(r,  cp)  by  the  relations  (2.5).  The  snail 
parameter  X,  which  characterizes  the  deviation  of  the  nonlinear  law  of  elastici¬ 
ty  from  Hooke’s  law,  has  the  dimension  1/bar2  and  a  magnitude  of  the  order  of 
10_s  to  10-6 ,  while  the  dimensionless  constant  gs ,  for  certain  nonferrous  metals 
and  their  alloys,  is  of  the  order  of  l(f  to  106  (Bibl.2). 

The  components  of  the  displacements  u(r,  cp)  and  v(r,  9)  for  the  nonlinear 
law  of  elasticity  (2.1)  under  the  conditions  (2.3)  are  determined  from  the  sys¬ 
tem  of  equations 
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^  =  ~9inc(F"+TF'  +  7rF")  +  w(“F"+2TF'+2^F") 
+  F  (2 M  +  w)  [n  ~TF'F”+7i {F'  -  F”F" + 3/?*> 

+  7^^  “  3  F.F„)  +  ~(fU  +  3F*)](-F„  +  2y  F,  +  2~F„) , 

tIt+t-  9^  (**+>*+>-)+ «**( 

+  ^-(2^  +  -^)  [*?r  -  ~  ft  ftr  +  ~  (fi  ~  ft,  F„  +  3ft„) 
+  2~(F,F„-3F,F„)  + 1(^  +  3fJ)j(2Frr  — i-fr 


(2.8) 


Thus,  the  solution  of  the  problem  reduces  to  an  integration  of  the  complex 
fourth-order  nonlinear  equation  (2.4)  or  of  the  system  of  nonlinear  equations 
(2.8)  under  definite  boundary  conditions  over  the  orifice  contour  as  well  as 
under  conditions  of  the  behavior  of  these  functions  ”at  infinity”  (Bibl.l). 


3.  Consider  orifices  of  a  shape  such  that  the  function 

Z*  =  F[C  +  e/(C)) ,  (z*  -  r'e* ;  r*  «  Rr ;  z  -  re* ;  f  =  Qe**) 


/1012 

(3.1) 


realizes  the  conformal  mapping  of  an  infinite  plane  with  a  circular  orifice  of 
unit  radius  onto  an  infinite  plane  with  an  orifice  having  the  shape  considered. 
In  the  function  (3.1) ,  R  is  the  true  constant  characterizing  the  dimensions  of 
the  orifice;  the  function  f(£)  depends  on  the  shape  of  the  orifice;  e  is  a  small 
parameter;  the  true  quantity  satisfying  the  condition  e  <  1  and  the  roots  of 
the  equation  1  +  ef*  (£)  =  0  should  lie  within  the  unit  circle  in  the  plane  Q. 

We  will  present  the  solution  of  eq.(2.4)  and  of  the  system  of  equations  (2.8)  in 
the  form  of  expansions  in  the  small  parameters  m-  and  e  (2.7): 


F(r,  <p;  A*;«)  =#o  ^  nkeJFtk’»(r,  <p) ; 


k-OJ-O 


u(r,r>  p\  e)  =  H0  £  /ikeJu(k-})(r,<p) , 

k-0j~O 
OO  OO 

v(r,  <p;  fi;  e)  =  H0  £  Mk^v(kJ)(r,  <p) , 

k-QJ-0 

where  Ho  is  selected  from  the  condition  H^i32/R4  =  1.  Hence, 


(3.2) 


(3.3) 
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(3.4.) 


The  stress  and  strain  components  in  the  coordinate  system  (p ,  9)  also  will  be 
presented  in  the  form  of  series  in  p  and  e: 


a»  =  2  2  HieJ<J(ek’J) ,  a,  -  nke‘otJ) , 

At-oy-o  /c-0y-0 

OO  CO 

r<x> =  2  2  A'1#"0; 

k-OJ-0 

u*  =  2  2  J) .  u»  -  2  £  fikeJu?>» , 


k-0j-0 


k^0j-0 


(3.5) 


(3.6) 


Substituting  the  function  F(r,  cp;  p;  e)  of  eq.(3*2)  into  the  fundamental  equa¬ 
tion  (2.4)  and  equating  to  zero  the  coefficients  in  the  presence  of  identical 
exponents  of  p* ,  eJ  ,  will  yield  the  equation  for  determining  the  function 
in  the  form  of 


AAF<k-J\r,  <p)  =  LkJ(F<0-0) . 


(3.7) 


We  will  present  explicit  expressions  for  the  right-hand  sides  of  eq.(3*7)>  for 
certain  values  of  k  and  j. 

For  k  =  0,  j=0,  1,  ... 

^oj(^<0,0)»  •••  =  0. 

For  k  =  1,  j  =  0 

£i(o(^M))  = 

The  developed  form  of  the  operator  Lo (r( 0,0 was  given  by  Tsurpal  (Bibl.B). 

For  k  =  2,  j  =  0 

L,, ,(!**•>,  ^  =  L1(F<0’,>, F<1,0)).  (3.10 

The  developed  form  of  the  operator  L^F^0’05  ,  F^1,0)  )  was  given  in  another  report 
(Bibl.9). 

For  k  =  1,  j  =1,  the  developed  form  of  the  operator  L^i  will  be 


(3.8) 

am 

(3.9) 
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where 


L,., (F^®),  F*0-1*)  =  7^>0,(-1t-  J*°;0)  +  |fi°’0)) 

+ ei)(^nr + ^ + -ino>i))  -  2[(y^,#) 

j_n«.«,)(i^o,i,  _  -  ^n#i  o))  (y  -  ^n°*i))] 


•  -|-zl  (I^0llMir(0'0)  +  r<O.0)J /««.!))  f 


(3.11) 


r(V°;»)  =  (/^r°'°))1  -  yjy-VjXO’O)  +  JL[(/v<°.°>)*  -  +  3(F<®-°>)*] 

+  ~  3F?®'®>F'®-®>)  +  +  3(F<®  •»>)*] ; 

J$;»{  =  2/*®'®>i*®-1>  -  +  jp-(o,i)/?(o.o)j  +  JL(2i*®.o>/*®.‘> 

-  /?(o.o)f’(o,n  _  /*n.i>/rfM>  +  efjo.ojfrt.D)  +  A  (/*«•»>  F^  + 

-  3F%-VP<?’»  -  3F$"F(W)  +  i(2F<y>F<y>  +  6F<0-0>F<°-1>) . 


(3.12) 


(3.13) 


Proceeding  analogously,  we  can  write  explicit  expressions  for  the  operators  Lkj 
at  any  (k  >  1,  j  >  l)  values  of  k  and  j.  The  solution  of  eq.(3»7)  is  sought  in 
the  form  of  a  Fourier  series 


=  2l{fi:)(r)oosm<p  +  g^J(r)  sinmy). 


(3.14) 


To  find  the  stress  components  Op ,  Og ,  and  TpQ  and  the  displacement  components 
Up  and  ue  in  the  curvilinear  orthogonal  coordinate  system*  (p ,  6 )  given  by  the 
function  (3*1),  we  utilize,  by  analogy  with  another  paper  (Bibl.6),  the  corre¬ 
sponding  formulas  for  conversion  from  the  polar  coordinates  (r,  cp)  to  the  curvi¬ 
linear  orthogonal  coordinate  system  (p,  9). 

Expanding  the  obtained  expressions  for  the  stress  components  arp,  <Jq  ,  and 
TpQ  and  for  the  displacement  components  Up  and  Uq  in  series  in  pi  and  e  and  /1014 
taking  into  account  the  form  of  the  function  (3*1),  we  have 


*  \  Q  5g  Q*  50* / 


*  The  linear  coordinate  p  =  1  coincides  with  the  contour  of  the  orifice  under 
study. 
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1  3 

1  3* 

W1 

e 

e *  30* 

L  J!!_\  _  nJ-m)  — L*|  0) ; 

>»  30*/  **  3g30  gj  h 


*  m-0  L 

+_L  Jl_  Jl]  j.  «/-•>  JliU 

g*  30*  3g*/  '  3g30  q  J 


+ 


MjLj 

1  3* 

3* 

(He1 

r  Q *  30* 

“3g* 

\F^k-m\Q,6; 


j- • 


W(M)  _  JjT  + 

»i“0 

y-i 

_  v(*.y)  JT1  [^y-")^*.")  _  ^y-")^*.")] . 


(3.15) 


(3.16) 


m—0 


Substituting  the  functions  F(r,  cp ;  M* ;  e)  [eq.(3.2)]  and  u(r,  cp) ,  v(r,  cp) 
[eq.(3.3)]  into  eq.(2.8)  and  equating  to  zero  the  coefficients  in  the  presence 
of  identical  exponents  of  pk  ,  eJ  ,  we  will  find  the  system  of  equations  for  de¬ 
termining  u(K,J)  (r,  cp),  v(k,J)  (r,  cp)  which  enter  into  eq.(3.l6). 

.  The  functions  F<k,J)  (p ,  9)  entering  into  eq.(3.15)  are  solutions  of  eqs.(3.7) 
F^  '  ‘  (r,  cp)  in  the  form  of  eq.(3.14),  in  which  the  variables  r  and  cp  are  re¬ 
placed  by  p  and  9 ,  respectively.  Then,  the  arbitrary  constants  entering  into 
fic,j  (r)  and  &,j  ,(n)  [eq.(3.14)]  are  determined  from  the  corresponding  boundary 
conditions  for  F;  k’  3  (p ,  9 ) ;  these  conditions  are  derived  from  the  expansions 
analogous  to  eqs.(3.2)  and  (3*3),  for  the  orifice  contour  under  study,  on  the 
basis  of  an  expansion  in  a  double  series  in  pi  and  e . 

The  stress  and  strain  coefficients  over  the  contour  of  the  orifice  in 
question  are  found  from  eqs.(3.4)>  (3*5),  and  (3.15)  at  p  =1. 

Consider  the  problem  in  which  the  stresses 


°t\r=  Vlr  “  Vi(.r,<p\p)  (3.17) 

are  specified  over  the  orifice  contour.  For  the  sake  of  universality,  let  us 
assume  that,  in  eq.(3*17)»  the  quantities  and  \|f3  depend  on  p .  The  equation 
of  the  orifice  contour  in  parametric  form  is  determined  by  the  function  (3.1)  /101$ 
and  may  be  written  as 

r  —  r(g,  0) ,  <p  =  <p(g,6)  ,  * 
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for  p  =  1. 


Having  utilized  eqs.(3*l)  and  (3.18),  we  will  present  the  right-hand  sides 
of  eq.(3«17)  as  double  series  in  p.  and  e 


\r  =  2  fike}y>{kJ)(0) ,  Te»|r  -  £  E •  (3.19) 

*-o  y-o  fc-oy-o 


Substituting  expressions  (3.4)  into  eq.(3.19),  assuming  that  p  =  1,  and  compar¬ 
ing  the  coefficients  with  identical  exponents,  we  obtain  the  relations 


<7)|r=V«0).  #'>| 


(3.20) 


From  eq.(3.15),  taking  eq.(3*20)  into  account,  we  obtain  the  boundary  conditions 
for  determining  the  (k,  jjth  function  F  }  (r,  cp)  in  the  form  of 


l~~ 

1  d 

i  a*  \ 

w 

~  q  dg 

e*  ae*/ 

d1  1 
dgdO  <j 


F«-»(q,0) 


(3.2i) 

-  -  4  *H..r  {£ K”  -  7 


+ 


F  w  -  £)}*'*  4-.‘ 


In  accordance  with  eq.(3.2),  the  solution  obtained  to  the  n^h  approximation  will 
be  construed  as  the  function 


k+J^n-l 

Fn(r,  r< =  ft 0 

k.J- 0 


HkeJF{kJ)(r,  <p) . 


(3.22) 


From  eqs.(3.15)  and  (3.21)  we  can  see  that,  for  each  function  F(lc’J)  ,  we  obtain 
the  boundary  problem  for  a  circular  orifice.  This  also  explains  why  all  the 
basic  equations  and  the  relations  (2.1),  (2.4),  (2.6),  and  (2.8)  were  written 
in  polar  coordinates. 

In  eqs.(3.15)  and  (3.16),  LiJ-l)  ,  ...,  LsJ-a)  are  differential  operators 
whose  form  depends  on  the  function  f(0  C eq. (3 «l)3 •  The  series  expressions  for 
these  operators  for  the  zeroth,  first,  and  second  approximations  are  given  by 
Savin  and  Guzf  (Bibl.6). 
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!+.  Consider  by  way  of  an  example  the  most  elementary  case,  namely,  /1016 
the  case  of  omnilateral  uniform  tension  produced  by  stresses  P  applied  to  an  in¬ 
finite  physically  nonlinear  isotropic  plate  with  an  elliptical  orifice  (Fig.l) 
obeying  the  elasticity  law  (2.1)  under  the  conditions  (2.3). 


The  function  (3 • l)  which  maps  the  exterior  of  the  elliptical  orifice  onto 
the  exterior  of  a  circle  of  unit  radius  for  this  case,  as  is  known  (Bibl.10), 
has  the  form 


where 


Z * 


-*[c+i]. 


K  -  «  = 


a~b 
a  +  b ' 


f  =  Q<t\  z*  =  Rre,<r 


(4-1) 


(4-2) 


and  where  a  and  b  are  the  semiaxes  of  the  ellipse  (Fig.l).  Obviously,  the 
f motion*  f(C)  in  eq.(3.l)  will,  in  this  case,  be  f(C)  =  1/C. 

The  approximate  solution  of  this  problem,  taking  into  account  the  three 
approximations,  reduces  to  the  successive  integration  of  eq.(3.7)  allowing  for 
the  form  of  the  operators  (3*8)  -  (3 .13). 

The  stress  functions  with  the  zero ^  F(0,0)  ,  first  F(1,0)  ,  and  second  F(2,0) 
approximations  for  the  mentioned  omnilaterally  stressed  physically  nonlinear 
elastic  plate,  with  a  circular  orifice,  are  known  (Bibl.9)  and  have  the  form 


pto.o)  —  _JL_(r*  _  2lnr) , 


(4.3) 


#  For  an  orifice  with  romded  angles,  the  fmction  f(C)  will  be  (Bibl.10),  for  a 
square  orifice,  f (C )  -  1/C3;  for  a  triangular  orifice,  f(C)  =  1/C2  ,  etc. 
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f(».o)  _ 


fit 


,o)  _  13,  ,  1  I 

ml  5  nr+  4  (' 


1  1  \  T 

/1017 

+  T7r)  +  lnrJ’ 

(4*4) 

1  1  47  1  51  1  \1 

3  r4  +  36  r*  +  80  r»'/  J 

(4.5) 

The  stress  function  for  a  linearly  elastic  plate  with  an  elliptic  orifice  has 
the  form  (Bibl.5): 


Substituting  the  function  (4*3)  into  eq.(3*12),  we  find  the  functions  T(0,0)  (r, 
cp)  in  the  form  of 

IW,-^(1+4)-  (1-7) 

,  Knowing  the  functions  (4*3)  and  (4*6)  we  find  from  eq.(3.13)  the  function 
T!oa)  (r,  ?): 


-  «££(l  -  s£  +  »£)«Or.  (4.e, 

Substituting  the  functions  (4*3),  (4*6),  (4*7),  (4-8)  and  their  derivatives 
into  eq.(3.1l),  we  find  the  series  expression  for  the  operator  L^i . 

The  differential  equation  (3.7)  for  the  function  F(1,1)  (r,  9)  will  be 

^(M,  + 64^(^  +  36^)0082,7  =  0.  (4.9) 

The  specific  integral  of  eq. (4.9)  becomes 

F~c'\r, <f)  =  -•£- (y  -V  +  y  — -)cos29? .  (4.10) 

The  general  integral  of  the  homogeneous  equation  (4*9)  will  be  taken  as 

OO 

Ftotyr,  (p)  =  (Cm3 r~m+2  +  Cmir-m)cosnup .  (4 . n) 
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The  integration  constants  C,3  and  0,4  in  eq.  (4. 11) are  determined  from  the 
boundary  conditions  (3 .21): 


(Ti+7^)F‘,',,foe,L,+ 6) 


=  0, 


e»i 


dqdO  q 


^■'Kq,  Q) 


=  0. 


(4.12) 


»“1 


Omitting  the  intermediary  calculations,  we  will  present  the  final  expressions 
for  the  function  F1’1*  (r,  cp) :  /1018 

Fihl\r,<p)  =  -±^(32-73i-+5-i-  +  36-l)cos2y.  (4-13) 


Let  us  consider  the  second  approximation  in  greater  detail.  The  stress 
function  Fs(r,  cp;  p,;  e)  [eq.(3*22)]  will  be 


Fj(r,  r,  /*;  e)  =  //o(f(0-0)  +  +  fF(0>1>  +  e’F(M)  +  /ifiF(lll)) . 


(4.14) 


From  eqs.(3«5),  taking  into  account  the  values  of  the  components  Opk,j  ,  aglc,J)  , 
and  TpQ,J)  of  eq.(3-15)  as  well  as  the  values  of  the  functions  F(  k,J)  (k  =  0,  1,  2, 
j  =  0,  1,  2)  in  eqs . (4*3)  -  (4-7),  (4*13),  and  (4.14),  we  will  determine  the 
stress  state  in  a  physically  nonlinear  thin  plate,  weakened  by  an  elliptical 
orifice,  to  a  second  approximation.  Over  the  orifice  contour,  the  coefficient 
of  stress  concentration  k(2)  =  ad/P  will  be* 


kw  =  =  2[1  -  1. 5002  p*  +  10. 6052V  +  2e  cos  20 

+  2e2  cos49  —  10,  6602c p1  cos20] . 


(4.15) 


5.  It  can  be  seen  from  eq.(4*15)  that,  if  allowance  is  made  for  the  physi¬ 
cal  nonlinearity  of  the  materials  satisfying  the  elasticity  relations  (2.1) 
under  the  conditions  (2.3),  the  coefficient  of  stress  concentration  will  non- 
linearly  depend  not  only  on  the  magnitude  of  the  tensile  stresses  P  (Fig.l)  and 
on  the  parameter  2.  (characterizing  the  strength  properties  of  the  plate  materi¬ 
al)  but  also  on  the  ellipticity  of  the  orifice,  characterized  by  the  parameter  e 
[eq. (4.2)].  Setting  e  =  0  in  eq.(4»15),  we  obtain  the  values  of  k  for  the  case 
of  a  circular  orifice  (Bibl.9).  Setting  2  =  0  in  eq.(4*15),  we  obtain  the 
values  of  k  found  by  Guz*  (Bibl.5)  when  using  the  mentioned  approximate  method 
for  the  case  of  an  elliptical  orifice  where  the  plate  material  obeys  HookeTs 
law.  For  this  last  case,  there  exists  an  exact  solution  (Bibl.10)  of  the  prob¬ 
lem.  A  comparison  of  the  corresponding  values  of  k  given  by  the  exact  solution 
(Bibl.10)  with  the  approximate  k  (Bibl.5)  will  yield  a  clear  idea  on  the  rate 


*  The  superscript  of  the  coefficient  of  concentration  in  eq.(4.15)  gives  the 
number  of  the  approximation. 
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TABLE  1 


nm 


a/b 

1.00 

1.05 

1.10 

1.20 

1.30 

130 

1.60 

Linear 

Theory 

Exact 

Solution 

2 

2 

2.101 

1.904 

2.212 

1.818 

2.444 

1.666 

2.616 

1.538' 

3  000 
1.333 

3.200 

1.250 

1-0 

App roximatc 
Solution 

2 

2 

2.097 

1.904 

2.198 

1.818 

2435 

1669 

2.587 

1.546 

2.960 

1360 

3.136 

1.289 

p 

In  the 

linear  theory , 

,  k  i» 

independent 

ofP 

Nonlinear 

Theory 

60 

1.920 

1.920 

2.002 

1.843 

2.084 

1.775 

2.248 

1.658 

2.412 

1365 

2.730 

1 .431" 

2.709 

1.557 

70 

1.904 

1.904 

1.980 

1.834 

2.055 

1.772 

2,208 

1.667 

2.361 

1384 

2.660 

1.469 

2.568 

1.667 

At 

80 

1.895 

1.895 

1.962 

1.832 

2.301 

1.777 

2.170 

1.686 

2.310 

1.616 

2387 

1.522 

2.414 

1.802 

90 

1.868 

1,868 

1.954 
1 .841 

2.014 

1-794 

2.138 

1.718 

2.264 

1.662’ 

2.517 

1393 

2.250 

1.966 

TABLE  2 


|  a/b 

1,00 

1,05 

1.10 

1.20 

1.30 

1.50 

1.60 

Linear 

Enact 

2.0000 

2.1010 

2.2120 

2.4440 

2.6160 

3.0000 

3.2006 

The ary 

Solution 

2.0000 

1.9047 

1.8182 

1.6667 

1.5387 

1.3333 

1.2509 

Approximated 

Solution 

2  0000 
2.0000 

2.0970 

1.9048 

2.1980 

1.8188 

2.4350 

1.66% 

23870 

1.5460 

2.9600 

1J600 

3.1360 

1.2896 

P 

In 

the  linear 

theory. 

k  is  independent  o 

f  k 

1.9606 

2.0406 

2.1308 

2  3102 

2.4874 

2.8298 

2.8987 

600 

1.9606 

1.8654 

1.7886 

1,6570 

1.5500 

1.3914 

1.4278 

1,9682 

2.0622 

2.1562 

23424 

2  5260 

2.8790 

2.9909 

1.9682 

1 ,8790 

1,7984 

1.6600 

1.5466 

1,3774 

1.3712 

1 .9226 

2,0050 

20878 

2.2534 

2.4196 

2.7380 

2.7248 

800 

1 ,9226' 

1.8450 

1.7756 

1 .6578 

1.5616 

1.4272 

13460 

1 .9474 

2.0368 

2.1262 

2-3044 

2.4914 

2.8200 

2.8810 

Theory 

1 .9474 

1 .8628 

1 .7868 

1.6564 

1.5504 

1 .3948 

1.4386 

1.9006 

1  9730 

2.0466 

2.1944 

2.3430 

2.6346 

2.5150 

1000 

1.9006' 

1 .8330 

1.7730 

1.6728 

1.5944 

1.4866 

1.7118 

1.9258 

2.0092 

2.0930 

2.2604 

2.4266 

2.7494 

2.7468 

, 

1.9258 

1,8472 

1.7768 

1.6572 

1.5612 

1.4222 

1.5302 

1200 

1 .8946 
1.8946 

1.9550 

1.8390 

2.0168 

1.7906 

2.1432 

1.7120 

2.2722 

1,6532 

2  3292 
13800 

2.2279 

1.9353 

1.9068 

1.9828 

2.0596 

2.2142 

2.3684 

2.6704 

23892 

*a 

1.9068 

1,8356 

1,7722 

1.6654 

13814 

1.4632 

1.6496 
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of  convergence  of  the  approximate  solution  of  the  problems  of  stress  concentra¬ 
tion  about  curvilinear  orifices  for  which  exact  solutions  are  lacking,  as  pro¬ 
posed  above.  Such  a  comparison  is  presented  in  the  first  two  rows  of  Tables  1 
and  2.  The  values  of  k  in  these  Tables  indicate  that  even  for  the  greatly  ex¬ 
tended  ellipse  a/b  =  1.6,  the  third  approximation  [eq.(4.15)3  gives  for  k  a 
very  good  agreement  with  the  exact  value  (the  difference  does  not  exceed 
2.5  -  3*0%).  Tables  1  and  2  present  values  of  k  [eq.(4.15)]  calculated  for  two 
points  A  and  B  on  the  orifice  contour  (see  Fig.l).  The  values  of  k  at  the 
point  A(3  =  0)  are  substituted  into  the  numerator  and  at  the  point  B(0  =  tt/2), 
into  the  denominator,  for  different  values  of  a/b,  P,  and  X.  The  values  of  X 

were  taken  from  another  report  (Bibl.7) J  Xx  =  1.02  x  10-5  — - -3-  (copper);  X3  = 

bar 

=  0.055  x  10~1 2 3 * 5 6  — 0" .  (aluminum  bronze);  X3  =  0.033  x  10~6  • — —  -5—  (open-  /1020 

bar  bar  - 

hearth  steel) .  The  numerical  data  presented  in  Tables  1  and  2  indicate  that : 
a)  the  ellipticity  of  the  orifice  -  as  in  the  classical  case,  i.e.,  when  the 
plate  material  obeys  Hooke’s  law  -  greatly  affects  the  coefficient  k  of  stress 
concentration;  b)  as  the  applied  tensile  stresses  P  (Fig.l)  are  increased,  k 
will  decrease  at  the  point  A  and  will  increase  at  the  point  B.  It  follows  that, 
as  the  numerical  values  of  the  parameters  P  and  X  increase,  a  consideration  of 
the  physical  nonlinearity  generally  yields  a  more  uniform  stress  distribution 
over  the  orifice  contour.  The  approximate  method  for  the  solution  of  the  prob¬ 
lems,  formulated  in  Point  1  above,  was  based  on  a  formal  expansion  of  the  re¬ 
quired  functions  in  double  series  over  the  small  parameters  p  and  e,  without 
evaluating  at  all  the  convergence  of  the  series.  An  idea  as  to  the  rapidity  of 
convergence  of  the  proposed  method  in  the  general  case  (Bibl.5j  9)  can  be  ob¬ 
tained  by  calculating  the  concentration  coefficients  k(  n-2)  ,  lr-  )k(n)  from 
formulas  analogous  to  eq.(4»15)  and  corresponding  to  the  stress  functions  F„-2> 
F„-i .,  Fn  for  the  preceding  approximations. 


BIBLIOGRAPHY 


1.  Savin,  G.N.:  Effect  of  Physical  Nonlinearity  of  Materials  on  Stress  Concen¬ 

tration  about  Orifices  (Vliyaniye  fizicheskoy  nelineynosti  materiala  na 
kontsentratsiyu  napryazheniy  okolo  otverstiy).  Prikl.  Mekh.,  Vol.9,  No.l, 

1963;  Vol.10,  No.l,  1964. 

2.  Kauderer,  G. ;  Nonlinear  Mechanics  (Nelineynaya  mekhanika).  Int.  Lit.,  1961. 

3.  Mors,  F.M.  and  Feshbakh,  G.;  Methods  of  Theoretical  Physics  (Metody  teoreti- 

cheskoy  fiziki) .  Int.  Lit.,  Vol.II,  1958. 

4*  Guz’,  A.N.;  Approximate  Solutions  of  Problems  of  the  Theory  of  Plates  and 
Slightly  Curved  Shells  for  Certain  Doubly-Connected  Regions  (0  pribli- 
zhennykh  resheniyakh  zadach  teorii  plastin  i  pologikh  obolochek  dlya 
nekotorykh  dvukhsvyaznykh  oblastey) .  Prikl.  Mekh.,  V0I.9,  No.l,  1963. 

5.  Guz’,  A.N.:  Approximate  Method  of  Determining  Stress  Concentration  about 

Curvilinear  Orifices  in  Shells  (0  priblizhennom  metode  opredeleniya 
kontsentratsii  napryazheniy  okolo  krivolineynykh  otverstiy  v  obolochkakh) . 
Prikl.  Mekh.,  Vol.8,  No. 6,  1962. 

6.  Savin,  G.N.  and  GuzT,  A.N.:  Stressed  State  about  Curvilinear  Orifices  in 


13 


Shells  (0  napryazhennom  sostoyanii  okolo  krivolineynykh  otverstiy  v 
obolochkakh) .  Izv.  Akad.  Nauk  SSSR,  OTN,  Mekh.  i  mashin.,  1964. 

7.  Tsurpal,  I.A.:  Stress  Concentration  about  a  Circular  Orifice  in  a  Physically 

Nonlinear  Elastic  Plate  in  the  Presence  of  Pure  Shear  (Kontsentratsiya 
napryazheniy  okolo  krugovogo  otverstiya  v  fizicheski  nelineynoy  uprugoy 
plastinke  pri  chistom  sdvige).  Prikl.  Mekh.,  Vol.8,  No. 4,  1962. 

8.  Tsurpal,  I. A.:  Approximate  Solution  of  the  Problem  of  the  Elastic  Equili¬ 

brium  of  a  Physically  Nonlinear  Elastic  Plate  with  a  Reinforced  Circular 
Orifice  (Priblizhennoye  resheniye  zadachi  ob  uprugom  ravnovesii  fizicheski 
nelineynoy  uprugoy  plastinki  s  podkreplennym  krugovym  otverstiyem) . 

Dokl.  Akad.  Nauk  Ukr  SSR,  No.l,  1963. 

9.  Tsurpal.  I.A.:  Approximate  Solution  of  Physically  Nonlinear  Plane  Problems 

of  Stress  Concentration  about  Orifices  (Priblizhennoye  resheniye  ploskikh 
fizicheski  nelineynykh  zadach  kontsentratsii  napryazheniy  okolo  otverstiy) . 
Prikl.  Mekh.,  V0I.9,  No. 6,  1963. 

10.  Savin,  G.N.:  Stress  Concentration  about  Orifices  (Kontsentratsiya  naprya¬ 
zheniy  okolo  otverstiy) .  GITTL,  1951  . 


14 


NASA-Langley,  1966 


F-408 


